Measurement of the strong coupling as in e + e -annihilation 
using the three-jet rate 



Jochen Schieck 3 '* 

"Ludwig-Maximilians-Universitdt Munchen, Am Coulombwall 1, D-85748 Garching, Germany 
and Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching, Germany. 



Abstract 
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We present a measurement of the strong coupling as using data collected with the JADE detector at centre-of-mass 
energies between 14 and 44 GeV. The three-jet rate as a function of the transition parameter y cut is determined using 
the Durham jet algorithm and the distribution is compared to QCD predictions. Recent theoretical calculations predict 
the three-jet rate at next-to-next-to-leading order. For the first time a measurement of as is presented using QCD 
predictions at next-to-next-to-leading order matched to predictions at next-to-leading logarithmic approximation, with 
subleading terms being included as well. We obtain 

a s (M z o) = 0.1 199 + 0.0010(stat.) ± 0.0021(exp.) ± 0.0054(had.) + 0.0007(theo.), 
being consistent with the world average value of as . 
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1. Introduction 

The strong coupling as is one of the fundamental 
parameters of the Standard Model of particle physics. 
Several possibilities exist to determine this parameter, 
one method being presented in this paper. In e + e~- 
annihilation with quarks in the final state hard gluons 
can be emitted from the quark with the emission prob- 
ability being proportional to the strong coupling as- 
Events with a hard gluon are identified by observing 
at least three particle jets in the detector. The the- 
ory of strong interaction, Quantum Chromo Dynamics 
(QCD), allows to predict the rate of three-jet events as 
a function of a single parameter, the strong coupling 
el's. Recent theoretical improvements allow to cal- 
culate the three-jet rate at next-to-next-to-leading-oder 
(NNLO) and a first measurement with data collected by 
the ALEPH detector at 91 GeV using the three-jet rate 
at a single resolution parameter was published [5]. 

In this paper we present the measurement of as using 
data collected with the JADE detector between the years 
1979 and 1986. For the first time NNLO predictions 
are matched with next-to-leading logarithmic approxi- 
mations (NLLA), including subleading soft logarithms, 
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to measure as with the three-jet rate. The complete de- 
scription of the analysis can be found in [6J]. 



The three-jet rate determined with the Durham 
clustering scheme 



In order to count the number of jets in the final state 
the reconstructed particle tracks and neutral clusters 
have to be clustered to jets. Several clustering schemes 
are available and we use the Durham clustering algo- 
rithm [7]. The number of jets depends on the reso- 
lution parameter y cut , a measure indicating the separa- 
tion between the second and the third reconstructed jet. 
The fixed order prediction of the three-jet rate deter- 
mined with the Durham clustering scheme is matched to 
NLLA predictions using the /^-matching scheme. Parts 
of the subleading soft logarithms, which are not in- 
cluded in the NLLA calculations, can be systematically 
controlled and included by adding a renormalisation 
scheme dependent K-term (8J. Including the K-teim 
significantly improves the description of the data. QCD 
predictions of the three-jet rate, without the expected 
contribution from hadronisation effects, are shown in 
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Figure 1 : QCD predictions for partons clustered to three jets using the Durham 
algorithm evaluated at a centre-of-mass energy of 35 GeV and the strong cou- 
pling being set to o , s(M z o)=0. 1 180. The predictions do not contain correc- 
tions expected from hadronisation effects. The figure shows the fixed order 
prediction (NNLO) together with the fixed order prediction matched to NLLA 
(NNLO+NLLA) and the matched predictions taking subleading soft logarithms 
into account (NNLO+NLLA+K). The band reflects the theoretical uncertainty 
estimated by setting the renormalisation scale factor Xp to 0.5 and 2. 



3. The JADE Detector, Data and Monte Carlo sam- 
ple 



The JADE detector collected data in e + e~- 
annihilation at the PETRA accelerator at DESY 
between 1979 and 1986. We use data taken at 14 GeV, 
22 GeV, 34.6 GeV, 35 GeV, 38.3 GeV and 43.8 GeV, 
adding up to a total integrated luminosity of about 195 
pb -1 . The number of selected hadronic events ranges 
between about 1000 events at 14 GeV and 20000 events 
at 35 GeV. The selection of hadronic events is based 
on the event multiplicity, momentum imbalance and 
on the total visible energy. The events are corrected 
for acceptance effects and the QCD predictions for 
the hadronisation of partons to hadrons. For these 
correction procedures Monte Carlo (MC) samples 
generated with PYTHIA, HERWIG or ARIADNE are 
used, with the generators being tuned to events taken 
with the OPAL detector at LEP at a centre-of-mass 
energy of 91 GeV. The expected contribution from 
e + e~ — > bb-events, as predicted by simulation, is sub- 
tracted from the observed three-jet rate. The excellent 
description of the three-jet event rate by simulation is 
visualised in Fig. [2] 




Fi gure 2: The three-jet rate as a function of the resolution parameter _v cu[ for 
data taken at a centre-of-mass energy of 35 GeV. The histogram shows the three- 
jet rate in comparison to simulation using PYTHIA.HERWIG and ARIADNE 
Monte Carlos models. The insert shows the difference to the model normalized 
to the combined statistical and systematic uncertainty. 

4. Measurement of the strong coupling as 

The measured and corrected three-jet rate together 
with the matched NNLO+NLLA+K QCD predictions 
is used to determine the strong coupling as- A x 2 - 
value, obtained from the difference between the mea- 
sured three-jet rate and the QCD predictions, is calcu- 
lated and minimised. The QCD predictions are applied 
with the renormalisation scale factor set to the natural 
choice Xfj-l. Events simulated with the PYTHIA event 
generator are used to correct the predicted three-jet rate 
for hadronisation effects. The fit range is determined by 
requiring the correction from hadronisation and detector 
effects to be small. In addition the leading log contribu- 
tion in the QCD prediction is required to be well below 
unity. The correlations between the different y cut bins 
are determined using MC events. The result obtained at 
a centre-of-mass energy of 35 GeV from the fit to the 
three-jet rate is shown together with the three-jet rate 
distribution in Fig. [3] 

4.1. Systematic Uncertainties 

Several sources of systematic uncertainties are con- 
sidered for the measurement of the strong coupling as. 
The fit is repeated and the difference between the vari- 
ation and the default analysis is taken as the symmetric 
systematic uncertainty. The different sources of system- 
atic uncertainties are arranged in three categories: 
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Figure 3: Result of the fit from data taken at a centre-of-mass energy of 
35 GeV. The insert reflects the difference between the QCD prediction calcu- 
lated with the strong coupling returned from the fit and the measured distribu- 
tion, normalized to the combined statistical and experimental uncertainty. 

• Experimental uncertainties: The fit is repeated 
with a somewhat modified event selection, differ- 
ent reconstruction software, different correction for 
detector effects or a slight variation of the fit range. 
The systematic uncertainties are evaluated apply- 
ing these variations and added in quadrature. The 
main contribution comes from using different de- 
tector correction procedures and using different re- 
construction software. 

• Hadronisation uncertainties: Different MC mod- 
els are applied to adjust the QCD predictions for 
hadronisation effects. The largest systematic un- 
certainty comes from evaluating the hadronisation 
correction with the HERWIG event generator in- 
stead of the PYTHIA event generator, as used in 
the default method. 

• Theoretical uncertainties: In the default fit the 
renormalisation scale parameter x^= p./ is set to 
the natural choise x^-l. As a systematic variation 
x M is set to 0.5 and 2 and the larger difference to 
the default fit is taken as the theoretical systematic 
uncertainty. 

4.2. Results 

The measurement of as is performed at each energy 
point separately and the final result is obtained by com- 
bining the results of as at the various centre-of-mass en- 
ergies. Due to large hadronisation corrections the result 



from the fit to the data taken at 14 GeV is not included 
in the combination. The combined value for as evolved 
to the energy M z o is a s (M z o) = 0.1199 + 0.0010(stat.)± 
0.0021(exp.) ± 0.0054(had.) ± 0.0007(theo.), consistent 
with the world average value of 0.1184 + 0.0007 U. 
The measurements of as at the various energy points is 
shown in Fig. |4] The evolution of as measured with the 
three-jet rate as a function of the centre-of-mass energy 
is consistent with the expectation from QCD. 
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Figure 4: Measurements of the strong coupling a$ determined with the three- 
jet rate using the Durham jet-finder algorithm and NNLO order predictions. For 
the measurements between 14 and 44 GeV the QCD predictions are matched 
with NLLA including soft logarithms. The energy points at 34.6 and 35 GeV 
are combined to a single entry. The line corresponds to the world average value 
with the corresponding uncertainty [^]. 

4.3. Study of the renormalisation scale dependence 
Besides fixing the renormalisation scale x M to the nat- 
ural choice the fit is repeated with x M and as being free 
in the fit. The renormalisation scale obtained by the 
fit is consistently below unity, but with large statisti- 
cal uncertainties and being consistent with one. The 
large uncertainty can be understood by looking in more 
detail to the variation of as with respect to x^. The 
change in as is rather flat, leading to small theoret- 
ical uncertainties and to large statistical uncertainties 
for the renormalisation scale x M . The results obtained 
from the various energy points are combined to a sin- 
gle value, leading to a s (M z o) = 0.1204 + 0.0009(stat.)± 
0.0021(exp.)±0.0059(had.)±0.0008(theo.), again being 
consistent with the world average value. 

4.4. Measurement of a§ using fixed order predictions 
The fit is repeated with fixed order predictions only 

and x M being set to one or leaving x M as free parame- 
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ter in the fit. For the fit with x^ being set to the natural 
choice the^- 2 /d.o.f. is considerably worse than the one 
obtained with matched QCD predictions (almost a fac- 
tor four in the ^f 2 /d.o.f. for the fit to the data collected 
at 35 GeV). This points to missing higher order terms in 
the QCD prediction. In addition the fit result also shows 
an increased sensitivity to the fit range. Using matched 
NNLO+NLLA+K QCD predictions the dependency on 
the fit range is almost negligible, while for NNLO pre- 
dictions the choice of the lower y cut value in the fit range 
has a significant impact on the result. 

Leaving as as well as x M free in the fit leads to a sig- 
nificantly improved description of the fit with^ 2 /d.o.f. 
values similar to the one obtained from the fit with 
matched QCD predictions. However, the value obtained 
for the scale x M is rather low (around 0.2) and within 
the statistical uncertainty inconsistent with being one. 
The uncertainty on the fitted scale parameter x M is about 
one order of magnitude smaller than the uncertainty on 
the scale parameter obtained with matched predictions. 
This indicates that the sensitivity of as with respect to 
the renormalisation scale x M is increased compared to 
the fit with matched QCD predictions. 

5. Conclusions 

We present the first measurement of as deter- 
mined with the Durham three-jet rate using matched 
NNLO+NLLA+K QCD predictions. For this measure- 
ment hadronic events from e + e~-annihilation taken at 
centre-of-mass energies between 14 and 44 GeV are 
used. The result a s (M z o) = 0.1199 ± 0.0010(stat.) + 
0.0021(exp.) ± 0.0054(had.) ± 0.0007(theo.) is consis- 
tent with the world average value. The largest uncer- 
tainty originates from modeling the transition from par- 
ton level, as predicted from the perturbative QCD calcu- 
lations, to the hadron level, as observed by the experi- 
ment. The increased precision of perturbative QCD pre- 
dictions, together with large hadronic e + e~ event sam- 
ples, will allow to scrutinize MC generators modeling 
the transition from partons to hadrons. A measurements 
of as using fixed order predictions only cannot describe 
the data satisfactorily. A fit leaving as as well as x M 
free returns a good description of the data, but the fitted 
value of x M is inconsistent with one. 

Fig. [5] summarizes the result obtained here to pre- 
vious analyses using jet rates and differential jet rates. 
The value obtained for as using a fit to the ^23 distri- 
bution returns the almost identical result for as as mea- 
sured with this analysis. The measurements obtained 
at higher centre-of-mass energies return smaller uncer- 
tainties, since the size of the hadronisation correction 



decreases with increasing centre-of-mass energy. The 
measurement of as using the four-jet rate and data taken 
between 14 and 44 GeV returns a smaller hadronisation 
uncertainty, leading to a decreased overall uncertainty. 
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Figure 5: Summary of as measurements using the three-jet rate, the four-jet 
rate and the differential V23 distribution applied to e + e~-data taken between 14 
and 91 Ge V. The first line corresponds to the measurement summarized in this 
paper, the second line to 10], the third line to 1 1 111 . the fourth line to [12] and 
the last line to [^]. The vertical band indicates the world average value of as [^]. 
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